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■ Abstract 

We show that a triple- Regge parametrization of inclusive single diffraction agrees 
with the data in the following two domains: (a) £ > 0.03 at all t, (b) \t\ > 1 GeV 2 
at all £. Since the triple- Regge parametrization fails when applied to the full £— t 
range of the total single-diffractive cross section, we conclude that damping occurs 
Q H ' only at low-£ and low-|i|. We give a ("toy") parametrization of the damping factor, 

D(£), valid at low-|t|, which describes the da^ 0,1 / dt data at the ISR and roughly 
accounts for the observed s-dependence of c^ a ' up to Tevatron energies. However, 
an effective damping factor calculated for the CDF fitted function for d 2 a l ^ al / d^dt 
at y^s = 1800 GeV and \t\ = 0.05 GeV 2 , suggests that, at fixed-£, damping increases 
^ . as s increases. 

We conjecture that, in the regions where the triple- Regge formalism describes the 
data and there is no evidence of damping, factorization is valid and the 'Pomeron- 
flux-factor may be universal. With the assumption that the observed damping is due 
to multi-"Pomeron exchange, our results imply that the recent UA8 demonstration 
that the effective 'Pomeron trajectory flattens for \t\ > 1 GeV 2 is evidence for the 
onset of the perturbative 2-gluon pomeron. Our damping results may also shed 
some light on the self-consistency of recent measurements of hard-diffractive jet 
production cross sections in the UA8, CDF and ZEUS experiments. 
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1 Introduction 



The inclusive (inelastic) production of beam-like particles, known as single diffraction, 
as in: 



and its analogous pp and ep interactions, presents one of the most interesting phenomena 
in strong interaction physics. An observed "rapidity gap" (absence of particles in a range 
of rapidity) between X and p/ in the final state signifies that the entire (color singlet) 
residual momentum of the proton, with beam momentum fraction, £ = 1— x p , participates 
in the interaction between it and the second beam particle. This effect is described in 
terms of the exchange of the Pomeron Regge trajectory flU, which embodies the idea of 
"factorization" . The momentum transfer, t, and the beam momentum fraction, x p , of the 
final state proton, pf, "tag" the corresponding parameters of the exchanged Pomeron. 

Since x p ~ 1 is observed to be the most likely beam momentum fraction of the final- 
state Pf, correspondingly the most likely value of the Pomeron's momentum fraction in 
the proton, £, is near zero. Nonetheless, at current collider energies the squared-invariant- 
mass of the system X in Eq. [|, s' = £s to good approximation, can be quite large. This 
fact led to a proposal to study hard scattering in such interactions, as a means of 
determining if the Pomeron possesses an observable partonic structure. The observation 
of the predicted hard scattering || £|, |[] supported the notion that the Pomeron behaves 
like a quasi-real object inside the proton with an effective Pomeron flux factor. An open 
question is to what extent such a flux factor is universal; for example, is it independent 
of beam particle or center-of-mass energy, or are there regions of phase space where 
factorization breaks down, due to interference with other more complex phenomena (e.g. 
multiple-Pomeron-exchange) ? 

One of the long-standing theoretical problems in high energy hadronic interactions 
has been the understanding of s-channel unitarization in Pomeron-exchange (diffractive) 
interactions. Empirically, one finds that the total diffractive cross section, cx^f ai , in Reac- 
tion [I] and in the corresponding pp interaction, initially rises from threshold and tends to 
level off or "flatten" at high energy ||, whereas the dominant triple-Pomeron[]?| descrip- 
tion of these processes (see below) continues to rise and soon exceeds the total pp cross 
section. There is no built-in mechanism in the pure triple-Pomeron process to account 
for the observed flattening of cx*^ aZ , and hence avoid the violation of unitarity. 

Figure 1 displays the problem [8-17]. The s-dependence of the total cross section for 
React. [I], <J t °d al is shownQfor Feynman-x p > 0.95 (or £ < 0.05) of the final state proton or 
antiproton (the domain where Pomeron-exchange is dominant). <J t s °J' al rises sharply from 
its threshold at 11.3 GeV beam momentum and gently levels off to ~9 mb at the highest 
Fermilab energy. The solid curve in Fig. [I] is the triple- Regge prediction discussed below. 
At high energies, it is in complete disagreement with the measured cross section. 

1 It is conventional to quote <T*2' a/ for £ m in < £, < 0.05, because the experimental acceptance usually 
depends weakly on x p in this region and because the integrated background from non- 'Pomeron exchange 
and other sources in this region is small enough to be neglected. 



p + Pi 



X + p f 



+ c.c. 
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In the continuing theoretical efforts to satisfy s-channel unitarity [18-22], the words, 
screening, shadowing, absorption and damping are all used[^3|] to describe effects due 
to multiple Pomeron exchange (two-Pomeron-exchange is also an important component 
in understanding pp elastic scattering [24] at low-|t|). These calculations have had vary- 
ing degrees of success. Goulianos has taken a more pragmatic approach to satisfying 
unitarity and suggested that the integral of the Pomeron flux factor in a proton should 
saturate at unity above yfs ~ 22 GeV. 

In the present Letter, we find that damping is confined to the low-£, low-|£| region. 
We continue the analysis of the UA8 Collaboration [25] and demonstrate that there are 
regions either at larger £ or at larger t, where the available data are well described by the 
triple-Regge formula and therefore require no damping. 

It is thus clear that the damping function depends on both £ and t. We attempt 
to determine the ^-dependence at low— \t\ of an "effective" multiplicative damping factor 
which could account for the discrepancies between data and solid curve in Fig. [l|. However, 
we call it a "toy damping factor" for several reasons. First, there are large gaps in the 
available data in Fig. [I] and some inconsistencies, therefore making it impossible to find 
a unique function. Secondly, the processes which give rise to the observed damping may 
imply a breakdown of factorization, in which case a simple universal damping factor may 
not exist at low-£ and low-|t|. Finally, there is some evidence that the nature of the 
^-dependence may itself depend on s at our highest energies. 



2 Triple— Regge phenomenology 

We briefly summarize the relevant formula. The Mueller-Regge expansion |7]] for the 
differential cross section of React. [I] is: 

^ = E G ijk {t) • ^W-M*) ■ (s')^- 1 (2) 

where aij(i) is the Regge trajectory for Reggeon i. The sum is taken over all possible 
exchanged Reggeons. The Gijf.(t) are products of the various Reggeon-proton and triple- 
Reggeon couplings and the signature factors. 

There are two dominant terms in Eq. ^] at small £, namely ijk = PPP and VV7Z, 
where the first term corresponds to the triple-Pomeron process, and the second corre- 
sponds to other non-leading, C=+ trajectories (e.g., f'2) in the Pomeron-proton interac- 
tion, The PPP term increases with increasing s', whereas the PPP. term decreases with 
increasing s' . 

Because the Pomeron is the highest-lying trajectory, when i — j — Pomeron, 1 — 2a 
is negative and the differential cross section increases sharply as £ — ► 0. This corresponds 
to the empirical observation that the most likely momentum fraction of the Pomeron in 
the proton, £, is near zero. Thus, the sharp rise in the triple-Regge prediction of cr*J aZ in 
Fig. [I] is due to the kinematic fact that the minimum value of £ decreases with increasing 

S as ^, m i n = S min / S 
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For fitting to data, Eq. |2| has been rewritten ||25|| as: 



d 2 cr sd 
d£dt 

where, 



[K F^tf e bt e- 2a{t) ] ■ M(sT + R (*?*], (3) 



the left-hand bracket is taken as the Pomeron flux factor, Fp/ P (t,£), and the right- 
hand bracket (together with the constant, ctq) is the Pomeron-proton total cross 



total 



section, 

The two terms in er^*"' correspond to the (V) afc (°) _1 factorQin the PPP and VV1Z 
terms in Eq. |2|. Thus, cr^* has a form similar to that of real particle cross 
sections ||2T 



The products Kao and Ka^R are, respectively, the values of G-p-p-p(t) and Gpppit) 
at t = 0. 

The Pomeron trajectory, a(t), has been shown [p5| to become relatively flat for 
\t\ > 1 GeV 2 (see next section); therefore a quadratic term is added to the stan- 
dard linear trajectory 0, a(t) = 1.10 + 0.25 t + a"t 2 . The non-zero value of b 
in e bt compensates for the presence of the quadratic component in the Pomeron 
trajectory^. 

|Fi(t)| 2 is the standard Donnachie-Landshoff ||M|1 form-factor .0 Since it has never 
been shown to describe React. [I] at large t, the e bt factor also serves as a possi- 
ble correction. Thus, the product, \Fi(t)\ 2 e bt , carries the t-dependence of the 
in Eq. |2| and is assumed to be the same for both G-p-pv an d Gpp-ji- Physically, 
this means that the Pomeron has the same flux factor in the proton, indepen- 
dent of whether the Pomeron-proton interaction proceeds via Pomeron-exchange 
or Reggeon-exchange. 



3 Where is triple— Regge applicable ? 

We already know from the information in Fig. |I| that the dominant contribution to 
the total cross section, namely the data with small-£ and small-|£|, are not described by 
the triple-Regge formalism; a damping of the Pomeron flux with increasing s is certainly 
required in this region. However, we see no reason to suppose that the same damping 
must apply to the entire £—t domainQ, as proposed by Goulianos0. However, this issue 

2 At very large s', rescattering effects may lead to a logarithmic dependence on s' as well as to other 
complications Q . 

3 If, as Donnachie and Landshoff[^7j have suggested, aip^ al depends on momentum transfer, that 
dependence would also be absorbed into the e bt factor. 

4p W _ 4mg-2.8t 1 

5/ Pomeron-exchange dominates out to £ ~ 0.05 and contributes significantly to £ ~ 0.1. 
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can be resolved by using available data to determine if there are regions in the £-i plane 
where the formalism does apply; that is, where damping is not required. 

The UA8 collaboration has recently reported|25| a (successful) simultaneous fit of Eq. ^| 
to their data on React, [l] at y^i = 630 GeV and the extensive data sample of the CHLM 
collaboration at the CERN-ISR with a/s = 23.5 and 30.5 GeV. They use the values, 
Qjjfc(O) — 1 = 0.10 and -0.32, obtained for the Pomeron and f /A2 trajectories, respectively, 



in fits to real-particle total cross sections|26|, |28|, [29]. 

The four free parameters, Kctq, a", b and R, are determined by fitting Eq. plus an 
empirical background function of the form, Ae *^ 1 , to the combined ISR-UA8 data set in 
the rangeQ, £ = 0.03-0.10. The fitted parameters are: 



Ka 
a" 
b 
R 



0.72 ±0.10 mb GeV" 2 
0.079 ±0.012 GeV" 4 
1.08 ± 0.20 GeV" 2 
4.0 ±0.6 



This value of a" was independently confirmed ||25|| in fits to the ^-dependence in the peak 
region with £ < 0.03. 

These triple-Regge results can be used to predict the total cross section of React. [|, 
<r*^ a ', be integrating Eq. |3] over the entire t-range, as well as for £ min < £ < 0.05. This 
yields the solid curve in Fig. [I] and illustrates the discrepancy]^] with the experimental 



a 



total 
sd 



Fig. |2] from Ref. [35] shows the ISR[[3(]] and UA8 data[^5j data in the restricted region, 
0.03 < £ < 0.04, where the small (~15%) non-Pomeron-exchange background can be 
ignored (the background is smaller than the size of the dots in the figure and about the 
same magnitude as the systematic uncertainty in absolute cross sections). The similarity 
of gg! at £ = 0.035 at both ISR and SppS energies reflects the fact that o t ^ t p al (s l ) has 

nearly the same value at both \fs' = 5 and 118 GeV.[] The term, (s')~ a32 , in Eq. |3] makes 
this possible. The solid curves in Fig. |2| are fits to these data without a background term 
and yield values of the 4 parameters which are in excellent agreement^ with those given 
above from the £ = 0.03 — 0.10 fit, thus lending credence to the relibility and stability of 
the fits. 

We have found a second region in the £-t plane, at small-£ but large-|t|, where the 
triple-Regge formalism also describes the ISR and UA8 data — with no additional free 
parameters. Fig. ||] shows the high momentum-transfer part of o- l ° d tal for £ m j n < £ < 0.05 
and the limited |t|-range, 1.0-2.0 GeV 2 , plotted vs. s for the ISR and UA8 data. The 
solid curve is the prediction of Eq. ^ using the above parameters. The dashed curve is 

6 For £ > 0.03, there are no concerns about experimental resolution causing "spill-over" from the large 
peak at £ ~ 0. 

7 This arises because, at fixed £ and t in Eq. ||, is proportional to Cpp a/ (£s). 
8 The same four parameters from this fit are (0.67±0.08), (0.078± 0.013), (0.88 ±0.19) and (5.0±0.6), 
respectively. 
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obtained by decreasing the 5-parameter by la from its central value. We see that, in 
contrast with the situation for <r*^ aZ , the triple-Regge formula, Eq. fj[ accounts for the 
observed s-dependence of the total single diffractive cross section in the high— \t\ range, 
1.0-2.0 GeV 2 . The different shapes of the curves in Figs, [l] and the solid curve in Fig. [3] 
are due to the t-dependence of the Pomeron trajectory. 

We have thus demonstrated that damping depends on both £ and t; it only exists in 
the small-£, smali-|t| region, and is not required by data away from that region — either 
at larger £ or at larger \t\. This could explain why CDF|17| reports abnormally large 
backgrounds in triple-Pomeron fits to their (low-£, low-|i|) data at y/s = 1800 GeV. For 
example, at £ = 0.035 (and \t\ = 0.05 GeV 2 ) where normally 15-20% background is found, 
their fitted formula corresponds to non-Pomeron-exchange backgrounds of 51%. Such a 
result can be expected if a ^-dependent Damping factor is required, but is left out of the 
fit; the fitted (large) background term compensates for the wrong ^-dependence in the 
triple-Regge equation without damping. Since CDF only reports the fitted functions, we 
compare our prediction of with the sum of their "signal" and "background" terms. 

This sum corresponds to the solid bands in Fig. f|| the curves are the vs. t predictions 
of Eq. [|at £ = 0.035. We see that, at 1800 GeV, the prediction agrees to within la of the 
CDF result; at 546 GeV, there is also good agreement at the lowest value, although 
their fitted t-dependences at the two energies are not self-consistent. 



4 Empirical determination of damping at small (£,£) 

The i-dependence of the disagreement between triple-Regge and the measured cross 
section is best seen by comparing the predictions with the experimental values of da l °J[ al /dt 
plotted vs. t. This is done in Figs. | for eight ISR energies^ and in Fig. || at the SppS- 
Collider[]. 

The dashed and dotted curves in Figs. [5] and ^| are the (undamped) triple-Regge 
predictions for da/dt, calculated by integrating Eq. |3]over the range, £ m j„ < £ < 0.05 (for 
the dotted curve, b is decreased by la from its central value). 

At the ISR energies, where the triple-Regge prediction only exceeds the data by 
about 10-15% (see Fig. |l|), the differences between dashed curves and data in Fig. [| 
are hardly noticable, because the dot sizes are roughly similar to the discrepancies. At 
y/s = 630 GeV, however, the same effect is larger and highly visible. At that energy, we 
see that there is a gradual transition from the low-|i| region which dominates cx*^ a ', and 
where the experimental a l °l al is smaller than the (undamped) triple-Regge prediction, 
to the higher— 1 1 1 region where the predictions agree with the data. This seems to be a 
smooth transition over the |t|-range, 0.5-1.0 GeV 2 . The situation at the lower ISR ener- 
gies in Figs. H is similar but less pronounced. We conclude that the discrepancies between 
predictions and data are confined to the low-|t| region. 

9 Some of these data were obtained with unequal energies for the two beams 

10 The UA4 points come from two independent runs, one at high-/? and one at low-/3 which allowed 
them to span most of the available i-range. 
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Since, as noted above, the calculated rise in cr t s °J; al with increasing s is due to the fact 
that, kinematically, £ m j n decreases with increasing s, it seems natural to introduce an 
empirical damping factor, D(£), in Eq. |3| which suppresses small ^-values; i.e., we strike 
at the "heart" of the problem. D(£) will be unity everywhere except at small-£. 

For -D(£), we have tried a "toy" damping function which decreases from unity for 
£ < 0.008, following a quadratic function as shown in Fig. |7|. The parameters of the 
quadratic function are choosen to reproduce the leveling-off of cr^f at ISR energies in 
Fig. |l|. To account for the Tevatron and SppS points, an additional (steep) fall-off is 
needed for £ < 0.0002. We arbitrarily use a cubic formQ. The dashed curves on Fig. [TJ 
shows how such a function accounts reasonably wellQ for cr*^ at high energies. 

The solid curves in Figs. |5| and ^| are calculated from Eq. |3] multiplied by the above 
damping factor. As expected from Fig. [I], the effect of damping is very small at ISR 
energies, but increases with s. At the energy of the SppS, however, the damping is about 
a factor of 3 in the low t region, which dominates a*° d * . There is good agreement with 
the damping predictions at low-|£| in both Figs. || and || 

While, as explained above, the parameters of the quadratic term are chosen to agree 
with the departure of the ISR cross sections from the triple-Regge prediction in Fig. p], 
there seems to be no reason a-priori why this formulation should be valid at higher en- 
ergies. To clarify this point, we assume that the formula CDF fitted [0 to their data 
is a sufficient description of and compare its ^-dependence at \t\ = 0.05 GeV 2 with 

that of Eq. §. The band in Fig. shows the ratio of the CDF at 1800 GeV to the 
triple-Regge prescription, which can be interpreted as an empirical damping factor. This 
decreases from near unity at £ = 0.03 to about 0.5 at £ = 0.01, but is insufficient to 
account for the factor of ~5 required by the 1800 GeV cross section in Fig. [I]; therefore 
an additional (rapid) decrease in the damping factor must occur at smaller £, analogous 
to our toy damping factor discussed above. 

The CDF function thus indicates that, at larger s, the onset of damping occurs at 
increasingly larger ^-values. However, this does not invalidate the solid (damped) curve 
calculated with the above D(£) and shown in Fig. |6], because the overall damping calcu- 
lation is not sensitive to the details of D(£) in the larger ^-region. 



5 Conclusions 

We first summarize the key points of this Letter: 

• We have a triple-Regge parametrization of inclusive single diffraction which agrees 
with the data in two domains of the £— t plane: (a) £ > 0.03 at all t, (b) \t\ > 1 GeV 2 
at all £ (Figs. 0, |3] and[|). Since the triple-Regge parametrization fails when applied 

n D(£) = 7000£ - (2.86 • 10 7 )£ 2 + (3.62 • 10 10 )£ 3 . D(f) has a total of 3 free parameters, since the 
quadratic and cubic have identical slopes and magnitudes at £ = 0.0002. These were chosen to be the 
parameters of the quadratic function and the slope of the cubic function at £ = 0. 

12 The bump between the ISR and SppS energies is due to the interplay between the 2-component a 1 ^ 1 
and the damping function at t; m in- 
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to the full £—t range of the total single diffractive cross section (£ m j n < £ < 0.05 and 
all t), we can conclude that damping occurs only at low-£ and low-|i|. 

• We have given a parametrization of the damping factor, D(£), valid for t < 0.5 GeV 2 , 
which describes all the low-|t| da t °i al / dt data at the ISR and roughly accounts for 
the observed s-dependence of cr*^' ai up to Tevatron energies (Figs. [TJ, ^], ^| and|7|). 

• An effective damping factor calculated for the CDF fitted function for at y/s = 
1800 GeV and \t\ = 0.05 GeV 2 , suggests that, at fixed-£, damping increases as s 
increases (Fig. |7J). 

These results raise a number of issues: We can conjecture that, in the regions where the 
triple-Regge formalism describes the data and there is no evidence of damping, factor- 
ization is valid and the Pomeron-flux-factor may be universal. A systematic program of 
testing universality of the flux factor in these regions should be carried out in pp, pp and 
ep interactions. 

Our damping results may shed some light on the measurements of the fK quantity 
from the cross sections for diffractive jet production in React. [I] and its analogue ep reac- 
tion (K is the normalization constant in Eq. [3] and / measures violation of the momentum 
sum rule, when / ^ 1). Assuming that the Pomeron has dominant gluonic structure [31 
there are three measurements of fK: 

UA8 § fK= 0.30 ±0.10 



CDF [|| fK= 0.11 ±0.02 



ZEUS El fK = 0.37 ±0.15 



The fact that the UA8 data is at large- \t\ where there is no damping, whereas the CDF 
data is at small- |t| where the damping factor in the region of the jets is of order 0.50, 
could account for the difference between the UA8 and CDF fK values. 

In addition, despite the large errors, it is interesting that the UA8 and ZEUS values for 
fK are consistent. This might be expected, if there is no damping in ep collisions at bigh- 
ts 2 . Of course, at low-Q 2 , where the photon exhibits hadronic properties, multi-Pomeron 
exchange, and hence damping, may result in smaller values of fK. 

In order to further study the possible s-dependence of the effective damping factor, it 
would be very useful to make detailed measurements of single diffraction in pp interactions 
at RHIC energies. This would fill in the large gap getween ISR and SppS-collider energies 
seen in Fig. [1]. 

We note that UA8EB] offers as possible explanations of their observed flattening of 



the Pomeron trajectory for \t\ > 1 GeV , either that it is an effect of multiple-Pomeron 



exchange, or that it is evidence for the onset of the perturbative 2-gluon pomeron p4|, p5 
In view of our observation that damping is not required in this t-region, it seems that 
the perturbative Pomeron, explanation is more likely. It may therefore be interesting to 
study the 7-Pomeron cross section from t min up through the \t\ > 1 GeV 2 region. 
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Figure 1: of pp or pp interactions (with £ < 0.05) vs. y/s demonstrating the 

flattening of the cross section with energy (a factor of two is included to account for both 
hemispheres). The insert is a blow-up of the ISR energy range. The upper curve is the 
Triple-Regge prediction described in the text; the dashed curve shows the consequence 
of multiplying it by the "toy" damping factor discussed in the text. The lowest energy 
points (open circles) are from bubble chamber experiments [8-12]; followed by those from 
the ISR (solid circles |L3| and triangles [0]), the SppS-Collider (solid square [15|, and 
the Tevatron (inverted triangles ||l7f). 
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Figure 2: Differential cross section, vs £, for 3 ISR measurements [30] and UA8[^| 
(single-arm cross sections). The curves correspond to the fit described in the text [25 
The points are averages of data in the £-range 0.03-0.04 (the non-Pomeron-exchange 
background in the data points is about the same magnitude as the diameter of the dots). 
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Figure 3: a 1 ^ 01 for £ < 0.05 and \t\ = 1.0-2.0 GeV 2 vs. y/s (a factor of two is included 
to account for both hemispheres). The solid curve is the same Triple-Regge prediction 
used in Fig. [I] (where it is integrated over all t); the dashed curve is the same, but with 
the "6" -parameter decreased by la from its central value. The lowest energy points 
(closed circles) are from the ISR ||13|| ; the highest energy point (solid square) is from the 
SppS-Collider f25[ 
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Figure 4: Bands are the CDF differential cross sections at £ = 0.035, calculated from 
their fitted functionsjlTj] (single-arm cross sections); the band widths are ±lcr errors on 
their amplitudes (as explained in the text, their "signal" and "background" are added 
together). The curves are from the same calculations used for the curves in Fig. ||. 
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Figure 5: da/dt vs. \t\ (single-arm cross sections) with £ < 0.05 at eight ISR 13, 14 



energies. The numbers shown in each plot are their s values (GeV 2 ). The solid curves 
are the integrals of Eq. ||] with damping included, using the parameters given in the text. 
The dashed and dotted curves are calculated without damping; the dashed curve uses the 
central value of b, while for the dotted curve, b is decreased by la from its central value. 
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Figure 6: Inclusive differential cross section for protons in React. [I] for x p > 0.95, 
measured in experiment UA8 with y/s = 630 GeV. and in experiment UA4[|I3|, |TB| with 
\/s = 546 GeV (single-arm cross sections; the integral is 4.7 ± 0.35 mb, or 9.4 ± 0.7 mb 
for er^f"'). The solid curve is the integral of Eq. § with damping included, as explained 
in the text. The dashed and dotted curves are calculated without damping; the dashed 
curve uses the central value of b, while for the dotted curve, b is decreased by la from its 
central value. 
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Figure 7: The damping function referred to in the text. In the £-range, 0.0002-0.008, 
the function shown is a quadratic, D(£) = 1 - 7500(0.008 - £) 2 - For f < 0.0002, the 
function is made to drop quickly to zero following a cubic function, as described in the 
text. The band is the ratio of the CDF at |t| = 0.05 GeV 2 and y/s = 1800 GeV to 
the triple-Regge prediction described in the text. 
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